OU-HEP/091015 
LPSC09155 



Beyond the Higgs boson at the Tevatron: 
detecting gluinos from Yukawa- unified SUSY 



Howard Baer 1 !, Sabine Kramlfl, Andre Lessa 1 ^ 
Sezen Sekmen 3 !^ and Heaya SummyQ 



L Dept. of Physics and Astronomy, University of Oklahoma, Norman, OK 73019, USA 
2 Laboratoire de Physique Subatomique et de Cosmologie, UJF Grenoble 1, CNRS/IN2P3, 
INPG, 53 Avenue des Martyrs, F-38026 Grenoble, France 
3 Dept. of Physics, Florida State University, Tallahassee, FL 32306 



Abstract 

Simple SUSY GUT models based on the gauge group SO{ 10) require t — b — r Yukawa 
coupling unification, in addition to gauge coupling and matter unification. The Yukawa 
coupling unification places strong constraints on the expected superparticle mass spec- 
trum, with scalar masses ~ 10 TeV while gluino masses are much lighter: in the 300-500 
GeV range. The very heavy squarks suppress negative interference in the qq — > gg cross 
section, leading to a large enhancement in production rates. The gluinos decay almost 
always via three-body modes into a pair of 6-quarks, so we expect at least four 6-jets per 
signal event. We investigate the capability of Fermilab Tevatron collider experiments to 
detect gluino pair production in Yukawa-unified SUSY. By requiring events with large 
missing Et and > 2 or 3 tagged 6-jets, we find a 5<r reach in excess of m g ~ 400 GeV for 
5 fbr 1 of data. This range in m g is much further than the conventional Tevatron SUSY 
reach, and should cut a significant swath through the most favored region of parameter 
space for Yukawa-unified SUSY models. 
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1 Introduction 



There is now an ongoing huge effort at the Fermilab Tevatron pp collider to extract a Standard 
Model Higgs boson signal from a daunting set of SM background processes. While such an 
effort is to be lauded- and if successful would complete the picture provided by the Standard 
Model (SM)- we note here that an even bigger prize may await in the form of the gluino 
of supersymmetric (SUSY) models [TJ. Current searches from CDF and DO collaborations 
have explored values of rrig up to ~ 300 GeV within the context of the minimal supergravity 
(mSUGRA or CMSSM) model [2l[3]. Here, we show that Tevatron experiments should- with 
current data sets- be able to expand their gluino search much further: into the 400 GeV 
regime, in Yukawa-unified SUSY, which is a model with arguably much higher motivation than 
mSUGRA [3]. Since Yukawa-unified SUSY favors a light gluino in the mass range 300—500 GeV, 
with the lower portion of this range giving the most impressive Yukawa coupling unification 
[5H9], such a search would explore a huge swath of the expected model parameter space. 

Supersymmetric grand unified theories (SUSY GUTs) based upon the gauge group 5*0(10) 
are extremely compelling [TO] . For one, they explain the ad-hoc anomaly cancellation within 
the SM and 577(5) theories. Further, they unify all matter of a single generation into the 16- 
dimensional spinor representation "0(16), provided one adds to the set of supermultiplets a SM 
gauge singlet superfield Nf (i = 1 — 3 is a generation index) containing a right-handed neutrinoQ 
Upon breaking of 5O(10), a superpotential term / 3 ^M N .N^N^ is induced which allows for 
a Majorana neutrino mass M^. which is necessary for implementing the see-saw mechanism 
for neutrino masses [ITJ. In addition, in the simplest 5O(10) theories where the MSSM Higgs 
doublets reside in a 10 of 5O(10), one expects t — b — r Yukawa coupling unification in addition 
to gauge coupling unification at scale Q = Mqut [T2],[T3] . In models with Yukawa coupling 
textures and family symmetries, one only expects Yukawa coupling unification for the third 
generation [14"] . 

In spite of these impressive successes, GUTs and also SUSY GUTs have been beset with a 
variety of problems, most of them arising from implementing GUT gauge symmetry breaking 
via large, unwieldy Higgs representations. Happily, in recent years physicists have learned 
that GUT theories- as formulated in spacetime dimensions greater than four- can use extra- 
dimension compactification to break the GUT symmetry instead [15J. This is much in the 
spirit of string theory, where anyway one must pass from a 10 or 11 dimensional theory to a 
4-d theory via some sort of compactification. 

Regarding Yukawa coupling unification in 5O(10), our calculation begins with stipulating 
the b and r running masses at scale Q = Mz (for two-loop running, we adopt the DR reg- 
ularization scheme) and the t-quark running mass at scale Q = m t . The Yukawa couplings 
are evolved to scale Q = M SU sy, where threshold corrections are implemented [IB] , as we pass 
from the SM effective theory to the Minimal Supersymmetric Standard Model (MSSM) effec- 
tive theory. From Msusy on to Mgut, Yukawa coupling evolution is performed using two- loop 
MSSM (or MSSM+RHN) RGEs. Thus, Yukawa coupling unification ends up depending on the 
complete SUSY mass spectrum via the t, b and r self-energy corrections. 

In this work, we adopt the Isajet 7.79 program for calculation of the SUSY mass spectrum 
and mixings [TT]. Isajet uses full two-loop RG running for all gauge and Yukawa couplings and 

1 Here, we adopt the superfield "hat" notation as presented in Ref. [T]. 
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soft SUSY breaking (SSB) terms. In running from Mqut down to M wea k, the RG-improved 
1-loop effective potential is minimized at an optimized scale choice Q = ^/m^m^, which 
accounts for leading two-loop terms. Once a tree-level SUSY/Higgs spectrum is calculated, 
the complete 1-loop corrections are calculated for all SUSY/Higgs particle masses. Since the 
SUSY spectrum is not known at the beginning of the calculation, an iterative approach must 
be implemented, which stops when an appropriate convergence criterion is satisfied. 

Yukawa coupling unification has been examined in a number of previous papers [SHHK^lCnil 
HBH20]. The parameter space to be considered is given by 

m 16 , m 10 , Ml, mi/2, A), tan/3, sign(fi) (1) 

along with the top quark mass, which we take to be m t = 172.6 GeV [2T]. Here, m 16 is the 
common mass of all matter scalars at Mqut, mio is the common Higgs soft mass at Mqut 
and M|, parameterizes either D-term splitting (DT) [5]|9][T9] or "just-so" Higgs-only soft mass 
splitting (HS) [5][6]|20]. The latter is given by Tn 2 H = m± =F 2Mp. As in the minimal 
supergravity (mSUGRA) model, mi/2 is a common GUT scale gaugino mass, Aq is a common 
GUT scale trilinear soft term, and the bilinear SSB term B has been traded for the weak scale 
value of tan /3 via the EWSB minimization conditions. The latter also determine the magnitude 
(but not the sign) of the superpotential Higgs mass term \i. 

What has been learned is that t — b — r Yukawa coupling unification does occur in the MSSM 
for fi > (as preferred by the (g — 2)^ anomaly), but only if certain conditions are satisfied. 

• tan f3 ~ 50. 

• The gaugino mass parameter mi/2 should be as small as possible. 

• The scalar mass parameter mie should be very heavy: in the range 8 — 20 TeV. 

• The SSB terms should be related as A\ = 2m\ G = 4mf 6 , with A < (we use SLHA [2"2~| 
conventions). This combination was found to yield a radiatively induced inverted scalar 
mass hierarchy (IMH) by Bagger et al. [23] for MSSM+right hand neutrino (RHN) models 
with Yukawa coupling unification. 

• EWSB can be reconciled with Yukawa unification only if the Higgs SSB masses are split 
at Mqut such that m 2 Hu < m 2 H H The HS prescription ends up working better than DT 
splitting [121120]. 

In the case where the above conditions are satisfied, Yukawa coupling unification to within a 
few percent can be achieved. The resulting sparticle mass spectrum has some notable features. 

• First and second generation matter scalars have masses of order mie ~ 8 — 20 TeV. 

• Third generation scalars, and /x are suppressed relative to m 16 by the IMH mechanism: 
they have masses on the 1 — 2 TeV scale. This reduces the amount of fine-tuning one 
might otherwise expect in such models. 

2 An exception is the case of highly split trilinears [23] . 
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• Gaugino masses are quite light, with m g ~ 300 — 500 GeV, m~o ~ 50 — 80 GeV and 
m~± ~ 100 - 160 GeV. 

Xi 

Since the lightest neutralino of S'O(IO) SUSY GUTs is nearly a pure bino state, it turns out 
that its relic density Q~oh 2 would be extremely high, of order 10 2 — 10 4 (unless it annihilates 
resonantly through the light Higgs [6], which is the case only in a very narrow strip of the 
parameter space). Such high values conflict with the WMAP observation [23], which gives 

Q CDM h 2 = P cdm/Pc = 0.1099 ±0.0124 (2a). (2) 

where h = 0.74 ± 0.03 is the scaled Hubble constant. Several solutions to the S'O(IO) SUSY 
GUT dark matter problem have been proposed in Refs. p3|l8"ll2"6~| . The arguably most attractive 
one is that the dark matter particle is in fact not the neutralino, but instead a mixture of 
axions a and thermally and non-thermally produced axinos a. Mixed axion/axino dark matter 
occurs in models where the MSSM is extended via the Peccei-Quinn (PQ) solution to the 
strong CP problem [27]. The PQ solution introduces a spin-0 axion field into the model; if the 
model is supersymmetric, then a spin-| axino is also required. The SO (10) SUSY GUT models 
with mixed axion/axino DM can [8] 1. yield the correct abundance of CDM in the universe 
(where a dominant axion abundance is most favorable), 2. avoid the gravitino/BBN problem via 
m(gravitino) ~ mie ~ 10 TeV and 3. have a compelling mechanism for generating the matter- 
antimatter asymmmetry of the universe via non-thermal leptogenesis [25] . A consequence of 
the mixed axion/axino CDM scenario with an axino as LSP is that WIMP search experiments 
will find null results, while a possible positive result might be found at relic axion search 
experiments |29j . 

A more direct consequence of the Yukawa-unified SUSY models is that the color-octet 
gluino particles are quite light, and possibly accessible to Fermilab Tevatron searches. Under 
the assumption of gaugino mass unification, the LEP2 chargino mass limit that m~± > 103.5 

X\ 

GeV normally implies that m g ~ 430 GeV, quite beyond the Tevatron reach. However, in 
Yukawa-unified SUSY, the trilinear soft breaking term is large: Aq ~ 10 — 20 TeV. Such a 
large trilinear term actually causes a large effect on gaugino mass evolution through two-loop 
RGE terms, as illustrated in Fig. [TJ Here the left frame shows the gaugino mass evolution 
for the mSUGRA model point with (mo, mi/2, A , tan/3, fi) = (500 GeV, 157 GeV, 0,10,+), 
which has m~± = 103.5 GeV, with m g = 426.1 GeV. In the right frame, we show the gaugino 
mass evolution for Point B of Table 2 of Ref. [8], but with a slightly lower m^j value. This 
point has the following GUT scale input parameters: m 16 = 10000 GeV, m 10 = 12053.5 GeV, 
M D = 3287.12 GeV, m 1/2 = 34 GeV, A = -19947.3 GeV, tan/3 = 50.398 and /i > 0. (tan/3 is 
input as a weak scale value.) In this case, the gaugino mass evolution is strongly affected by the 
large two-loop terms, resulting in a much smaller splitting between gaugino masses M 2 and M3. 
Here, we find (after computing physical masses including one-loop sparticle mass corrections) 
that while m~± = 108.2 GeV, the gluino mass is only m g = 322.8 GeV. This value of m g may 
well be within range of Tevatron discovery, even while respecting chargino mass bounds from 
LEP2. 

In Yukawa-unified models, the b and r Yukawa couplings are large, while the top and 
bottom squark masses are much lighter than their first/second generation counterparts. As 
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Figure 1: Evolution of soft SUSY breaking gaugino mass parameters using two-loop RG evolu- 
tion in the case of mSUGRA and in the case of the HS model with parameters as listed in the 
text. 

a consequence, gluino decays to third generation particles- in particular decays to 6 quarks- 
are enhanced. In addition, gluino pair production via qq fusion is normally suppressed by fr- 
aud s-channel interferences in the production cross section. For trig ~ 10 TeV, the negative 
interference is suppressed, leading to greatly enhanced gluino pair cross sections. Use may be 
made of the large gluino pair production cross section, and the fact that each gg production 
event is expected to have four or more identifiable 6-jets, along with large E™ 1SS , to reject SM 
backgrounds. 

In this letter, we examine gluino pair production at the Fermilab Tevatron collider. While 
negative searches for gluino pair production have been made, and currently require (under 
an analysis with ~ 2 fb _1 of integrated luminosity) nig ~ 308 GeV [21 [3] in mSUGRA-like 
models, use has not yet been made of the large gluino pair production cross section and high 
6-jet multiplicity expected from Yukawa unified modelsjj Here, we point out the importance of 
exploiting the 6-jet multiplicity to maximize the reach. By requiring Tevatron events with > 
4 jets plus large E™ 1SS , along with > 2 or 3 tagged 6-jets, QCD and electroweak backgrounds 
can be substantially reduced relative to expected signal rates. We find that the CDF and DO 
experiments should be sensitive to m g ~ 400—440 GeV with 5 — 10 fb~ of integrated luminosity. 
Thus, Tevatron experiments are sensitive to much higher values of gluino mass than otherwise 
expected from conventional searches. With 5 — 10 fb" 1 of data, Tevatron experiments can 
indeed begin to explore a large swath of Yukawa-unified SUSY model parameter space. 

In Sec. [2X we review gluino pair production total cross sections and expected branching 
fractions, and introduce a special Yukawa-unified SUSY model line. In Sec. [31 we provide 
details of our event simulation program, and show how the requirement of events with > 4 jets 
plus large E^ 1SS , along with > 2 or 3 identified 6-jets, rejects much SM background, at little 

3 The utility of &-jet tagging for extracting SUSY signals at the LHC has been examined in Ref's [50] . 
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cost to signal. We provide our reach results versus rrig. In Sec. SJ we present a summary and 
conclusions. 

2 Production and decay of gluinos at the Tevatron 
2.1 Gluino pair production 

Recent studies of squark and gluino pair production at the Fermilab Tevatron collider, using 
data corresponding to 2 fb -1 of integrated luminosity and a beam energy of yfs = 1.96 TeV, 
have produced limits at the 95% CL that m g > 280 GeV (in the case of CDF [2]), and m g > 308 
GeV (in the case of DO [3]). These studies- in the parts focused on gluino pair production- 
essentially asked for the presence of events with > 4 hard jets, plus large E™ 1SS and large Ht, 
where Ht is the scalar sum of the Ets of all identified jets in the event, beyond an expected 
SM background level. These studies do not use some of the unique characteristics common to 
gluino pair production in Yukawa-unified SUSY, so we expect Tevatron experiments to be able 
to do much better in this case. 

First, we present the expected total cross section rates for gluino pair production in Fig. |2l 
displaying leading order (LO) and next-to-leading order (NLO) cross sections as given by 
Prospino [31] . We adopt a common first /second generation squark mass of rriq = 10 TeV, 
and take the Tevatron energy as y/s = 1.96 TeV. We see from the figure that for m g = 300 
GeV, the cross section is about 900 fb, dropping to about 65 fb for m g ~ 400 GeV. Moreover, 
it remains at the level of several fb even for m g as high as 500 GeV. 

These cross sections are well in excess of those which enter the CDF and DO search for 
gluino pair production. To understand why, we first note that gluino pair production for 
m g ~ 300 — 500 GeV is dominated by valence quark annihilation via qq fusion at the Tevatron. 
The gg fusion subprocess is dominant at much lower gluino masses, where the gluon PDFs 
have their peak magnitude at small parton fractional momentum x. The qq —> gg subprocess 
cross section receives contributions from s-channel gluon exchange, along with t- and u- channel 
squark exchange diagrams [32]. The st- and sw-channel interference terms contribute negatively 
to the total production cross section, thereby leading to an actual suppression of a(pp — > ggX) 
for rriq ~ rrig. For rriq ^> m g on the other hand, the t-channel, w-channel and interference terms 
are all highly suppressed, leaving the s-channel gluon exchange contribution unsuppressed and 
dominant. The situation is illustrated in Fig. [3j where we plot the LO and NLO gluino pair 
production cross section for m g = 300, 400 and 500 GeV versus rriq. We see that as rrig grows, 
the total production cross section increases, and by a large factor: for m g = 400 GeV, as rriq 
varies from 400 GeV to 10 TeV, we see a factor of ~ 10 increase in total rate! 

At the present time- Fall 2009- CDF and DO have amassed over 5 fb -1 of integrated 
luminosityjj Thus, if m g ~ 400 GeV, there could be ~ 300 gluino pair events in each group's 
data. Such a large event sample may well be visible if appropriate background rejection cuts 
can be found. The exact collider signatures depend on the dominant gluino decay modes, which 
we discuss in the next section. 

4 It is expected that CDF and DO will reach the ~ 10 fb^ 1 level during 2010. 
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Figure 2: Production cross section a(pp — > ggX) in fb at the ^/i = 1.96 TeV Fermilab Tevatron 
collider versus m- g , for rrig = 10 TeV. 




Figure 3: Cross section of gluino pair production (in fb) at the Fermilab Tevatron collider as a 
function of nig, for rrig = 300, 400 and 500 GeV. Dashed is LO QCD, while solid is NLO, as 
given by Prospino. 
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2.2 Gluino decays in Yukawa- unified SUSY 

To examine the gluino decay modes in Yukawa-unified SUSY, we will adopt a model line 
which allows us to generate typical Yukawa-unified models over the entire range of m g which 
is expected. First, we note in passing that Yukawa unification is not possible in the mSUGRA 
model, since the large t — b — r Yukawa couplings tend to drive the m 2 H soft SUSY breaking 
term more negative than m 2 H , in contradiction to what is needed for an appropriate breakdown 
of electroweak symmetry. Yukawa-unified models can be found if one instead moves to models 
with non- universal Higgs masses, where m 2 Hu < m 2 Hd already at the GUT scale [18j[33]. In this 
case, m 2 Hu gets a head start in its running towards negative values. Detailed scans over the 
parameter space in Ref. [6] using the parameter space in [1] found a variety of solutions in the 
Higgs splitting (HS) model. We will adopt Point B of Table 2 of Ref. [8] as a template model. 
This point has the following GUT scale input parameters: m 16 = 10000 GeV, m 10 = 12053.5 
GeV, M D = 3287.12 GeV, m 1/2 = 43.9442 GeV, A = -19947.3 GeV, tan/3 = 50.398 and 
yU > 0, (where tan/3 is again at the weak scale). The Yukawa couplings at Mqut are found to 
be ft = 0.557, fb = 0.557 and f T = 0.571, so unification is good at the 2% level. The gluino 
mass which is generated is rrig = 351 GeV. 

If we now allow m 1( / 2 to vary, we still maintain valid Yukawa- unified solutions over the range 
of mi/2 : 35 — 100 GeV, corresponding to a variation in m g : 325 — 508 GeV. (The Yukawa 
unification gets worse as m\ji increases, and at m\ji = 100 GeV diminishes to 7.3%.) The 
value of the chargino mass at m\/2 = 35 GeV is m~± = 108 GeV, i.e. slightly above the LEP2 
limit. We will label Point B with variable as the Higgs splitting, or HS, model line. The 
value of the light Higgs boson is ~ 127 GeV all along the HS model line. 

Armed with a Yukawa-unified SUSY model line, we can now examine how the gluino decays 
as a function of gluino mass. The gluino decay branching fractions as calculated by Isajet are 
shown in Fig. HI Here, we see that at low rrig ~ 325 GeV, the mode g — > bbx® occurs at over 
60%, and dominates the g — > bbxi branching fraction, which occurs at typically 10-20% [34] . 
As rrig increases, the decay modes g — > tbxl + c - c - grows from the kinematically suppressed 
value of below 10% at m g ~ 325 GeV, to ~ 40% at m g ~ 500 GeV. All these dominant decay 
modes lead to two 6s per gluino in the final state, so that for gluino pair production at the 
Tevatron, we expect collider events containing almost always > 4 jets +E^ 1SS , with > 4 6-jets. 
Even more fr-jets can come from \o. decays, since X2 ~^ bt>Xi a ^ around 20% all across the HS 
model line. Only a very small fraction of gluino decays, less than 10%, lead to first/second 
generation quarks in the final state. 

3 Reach of the Fermilab Tevatron for gluinos in Yukawa- 
unified SUSY 

Next, we examine whether experiments at the Fermilab Tevatron can detect gluino pair pro- 
duction in the HS model line assuming 5-10 fb" 1 of integrated luminosity. We generate signal 
and background events using Isajet 7.79, with a toy detector simulation containing hadronic 
calorimetry ranging out to \r/\ < 4, with cell size Ar/ x A0 = 0.1 x 0.262. We adopt hadronic 
smearing of AE = 0.7/VE and EM smearing of AE = 0.15/VE. We adopt the Isajet GET- 



7 



350 400 450 500 

m_ (GeV) 



Figure 4: Dominant gluino branching fractions versus m g for the Yukawa-unified SUSY HS 
model line. 

JET jet finding algorithm, requiring jets in a cone size of AR = 0.5 with > 15 GeV. 
Jets are ordered from highest Et (ji) to lowest Et- Leptons within \r)t\ < 2.5 (£ = e, /i) are 
classified as isolated if pt(£) > 10 GeV and a cone of AR = 0.4 about the lepton direction 
contains Et < 2 GeV. Finally, if a jet with \rjj\ < 2 has a I?-hadron with Et > 15 GeV within 
AR < 0.5, it is tagged as a 6-jet with an efficiency of 50%. Ordinary QCD jets are mis-tagged 
as 6-jets at a 0.4% rate [35] . 

We also generate SM background (BG) event samples from W + jets production, Z + bb pro- 
duction, ti production, vector boson pair production, hadronic bb production, bbbb production, 
tibb production and Zbbbb (followed by Z — > vv) production The W + jets sample uses QCD 
matrix elements for the primary parton emission, while subsequent emissions (including g — >• bb 
splitting) are generated from the parton shower. For Z + bb, we use the exact 2 — > 3 matrix 
element, which is pre-programmed into Isajet using Madgraph [36]. We use AlpGen [37] plus 
Pythia [38J for bbbb and tibb production, and Madgraph plus Pythia for Zbbbb production [36J . 

For our first results, we exhibit the distribution in £"™ 1SS in Fig. \5\ as generated for the HS 
model line Pt. B (with m g = 350 GeV) as the blue histogram, along with the summed SM 
backgrounds (gray histogram). While the signal histogram is harder than the BG histogram, 
the BG level is so high that signal doesn't exceed BG until E^ ms ~ 300 GeV. Of course, this 
Pt. B gluino mass is well beyond the current Tevatron gluino mass limits, so this is easy to 
understand. 

To do better, we must adopt a set of cuts that selects out canonical gluino pair production 
events. Here, we will follow the recent papers BMPT [39], CDF [2] and DO [3], and require the 

5 We do not take into account the QCD dijet backgrounds which turn out to be negligible after the cuts 
described below. 
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Figure 5: Distribution in missing Et from gluino pair production in Yukawa-unified SUSY 
Pt. B, along with summed SM backgrounds (gray histogram), with minimal cuts listed in text. 

cuts listed in Table [TJ 



cuts 


£miss 


Hj< 


E T {jl) 


E T (j2) 


E T (jS) 


E T {jA) 


BMPT 


> 75 GeV 




15 


15 


15 


15 


CDF 


> 90 GeV 


280 


95 


55 


55 


25 


DO 


> 100 GeV 


400 


35 


35 


35 


20 



Table 1: Sets of cuts from Ref's 13^ , JM/ and /3]/ used in this analysis. In addition we require 
throughout > 4 jets, no isolated leptons, at least one jet with \r]j\ < 0.8 and A0(jT,j'2) < 160°. 

We have yet to make use of the high 6-jet multiplicity which is expected from Yukawa- 
unified SUSY. In Fig. El we plot the multiplicity of 6-jets expected from SM background 
(brown histogram), and the summed BG plus signal from HS Pt. B. (The BG in the = 
channel is very under-estimated, since we leave off QCD multi-jet production.) We see that the 
BG distribution has a sharp drop-off as increases. Especially, there is a very sharp drop-off 
in BG in going from the = 2 to the = 3 bin. When we add in the signal distribution, we 
see the histogram expanding out to large values of rib due to the presence of 4-6 6-jets per SUSY 
event. For the softer BMPT cuts, the signal hardly influences the n& = 0, 1, 2 bins. However, 
in the rib = 3 bin, there is a huge jump in rate, reflecting the presence of a strong source of > 3 
6-jet events. In the case of the CDF and DO cuts, which are much harder, the total BG is much 
diminished. In this case, the summed signal plus BG distribution actually becomes rounded, 
and is again much harder than just BG alone. For the CDF (DO) cuts, signal exceeds BG in 
the ri;, = 2 bin by a factor of 2 (3). By the time we move to the = 3 bin, then for both 
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Figure 6: Distribution in tagged b-jet multiplicity for gluino pair production in Yukawa-unified 
SUSY Pt. B, along with summed SM backgrounds (gray histogram), after cut sets BMPT, CDF, 
and DO given in Tabled . 

CDF and DO cuts, signal exceeds BG by over an order of magnitude. Using soft cuts and low 
6-jet multiplicity, one should gain a good normalization of total BG rates. Then, as one moves 
towards large 6-jet multiplicity nf, > 2 or 3, there should be much higher rates than expected 
from SM BG. 

Table [2] shows a listing of expected contamination from each BG source after the different 
sets of cuts. The hard E^ 1SS and Ht cuts largely eliminate the bb BG. The isolated lepton veto 
combined with large E^ 1SS requirement cuts much of W + jets. The remaining large BGs come 
from tt production and Z + bb production, where Z — > vv. Requiring > 4 jets along with large 
Ht for the CDF and DO cuts largely reduces Zbb to small levels, leaving tt as the dominant 
BG. 

In light of these results, we proceed by requiring BMPT, CDF or DO cuts, along with 
• rife > 2 or 3. 

In Fig. [7] we plot the resultant SM background (blue dashed lines), along with expected signal 
rates for the HS model line (full lines) for the three sets of cuts with > 2 (upper row) as 
well as Uf, > 3 (lower row). The SM background comes almost entirely from tt production. 
The third 6-jet in ti production can come from additional g — >■ bb radiation, or from QCD 
jet mis-tags. Since the dominant BG comes from tt production, and the a(pp — > tiX) cross 
section is well-known from standard top search channels, the background should be rather well 
understood. 

We see from Fig. [7J that signal actually exceeds BG for a substantial range of m~ g for all 
cases except the BMPT cuts with n& > 2. We also compute the signal cross sections required 
for a 5a discovery for each selection assuming 5 and 10 fb^ 1 of integrated luminosity, shown as 
the dot-dashed and dotted lines, respectively. The significance in as is derived from the p- value 
corresponding to the number of S+B events in a Poisson distribution with a mean that equals 
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BG 


<r(fb) 


events 


BMPT 


CDF 




DO 








> 26s 


> 36s 


> 2bs 


> 36s 


> 26s 


> 36s 


66 


3.8 x 10 8 


10 6 


— 


— 


— 


— 


— 


— 


ti 


5.9 x 10 3 


10 6 


51.9 


1.3 


8.6 


0.3 


3.9 


0.14 


bb+{Z -»■ vv) 


1.3 x 10 4 


10 6 


15.7 


0.4 










W + jets 


4.8 x 10 6 


5 x 10 6 


1.9 












VV 


9.7 x 10 3 


10 6 


0.6 


0.01 










bbbb 


6.3 x 10 4 


9.7 x 10 5 


0.39 


0.13 


0.065 


0.065 






ttbb 


11 


4.1 x 10 5 


0.39 


0.13 


0.066 


0.019 


0.037 


0.013 


bbbb + (Z z/z/) 


0.54 


6.6 x 10 3 


0.03 


0.01 


< io- 2 


< 10- 2 


< 10~ : 


1 < 10" 3 


Total 


70.9 


1.98 


8.7 


0.38 


3.94 


0.15 



Table 2: SM backgrounds in fb before and after cuts BMPT, CDF and DO for rib > 2 and > 3. 
The pt range for bb subprocess generation is 15 — 200 GeV. The px range for ti subprocess 
generation is 10 — 300 GeV. The yl range for Zbb subprocess generation is 100 — 400 GeV. In 
the above, V = W or Z. 



to the number of background events. The best reach is achieved with the hard DO cuts. In this 
case, requiring nb > 2, we find that signal exceeds the 5a level for 5 (10) fb -1 of integrated 
luminosity for m g < 395 (410) GeV. Requiring n& > 3, the 5a reach for 5 (10) fb -1 increases to 
m g = 405 (430) GeV. Thus, in the case of Yukawa-unified SUSY where an abundance of 6-jets 
are expected to accompany gluino pair production, we expect Fermilab Tevatron experiments 
to be able to probe values of m g to much higher values than have previously been found. 

Since the value of m g is expected to lie in the range 300-500 GeV for Yukawa-unified 
models, and in fact the Yukawa unification is best on the lower range of m g values, it appears 
to us that CDF and DO, using current data samples, stand a good chance of either discovering 
Yukawa-unified SUSY, or excluding a huge portion of the allowed parameter space. 

4 Conclusions 

In this paper, we explored the capability of the CDF and DO experiments to search for gluinos 
with properties as predicted by supersymmetic models with t—b—r Yukawa coupling unification. 
While a vast effort is rightfully being placed by CDF and DO to search for the SM Higgs boson, 
a potentially bigger prize- the gluinos from supersymmetric models- could be lurking in their 
data. The Yukawa-unified SUSY model is extremely compelling, in part because it combines 
four of the most profound ideas in physics beyond the SM: 5O(10) grand unification (which 
unifies matter as well as gauge couplings), weak scale supersymmetry, see-saw neutrino masses 
and the Peccei-Quinn- Weinberg- Wilczek solution to the strong CP problem. While we do not 
present a specific model which incorporates all these ideas into a single framework, a wide array 
of low energy, collider and astrophysical data give some indirect and also direct support to each 
of these ideas. The requirement of Yukawa coupling unification forces upon us a very specific 
and compelling sparticle mass spectrum, including first/second generation scalars at the ~ 10 
TeV scale, while gluinos are quite light: in the ~ 300-500 GeV range. We investigated here 
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Figure 7: Reach of the Fermilab Tevatron collider for gluino pair production in Yukawa-unified 
SUSY HS model line, versus m g . We show the reach for 5 and 10 fb~ x of integrated luminosity. 



whether these light gluinos are accessible to Tevatron searches for supersymmetry. 

Our main result is that the CDF and DO experiments should be already sensitive to gluino 
masses far beyond currently published bounds (which lie around the 300 GeV scale). This is 
due to three main factors: 

1. Two-loop RGE effects allow for gluinos as light as 320 GeV in the Yukawa-unified model 
with multi-TeV trilinear soft terms, even while respecting LEP2 limits on the chargino 
mass. In the case of generic SUSY models with TeV scale soft parameters, the LEP2 
chargino mass limit usually implies m g ~ 425 GeV. 

2. Gluino pair production cross sections with m g ~ 300 — 500 GeV are enhanced at the 
Tevatron due to the extremely high squark masses expected in Yukawa-unified SUSY. The 
huge value of nig acts to suppress negative interference effects in the qq — > gg subprocess 
cross section, leading to elevated production rates. 

3. Gluinos of Yukawa-unified SUSY decay through cascade decays to final states almost 
always containing four 6-jets, and sometimes six or eight 6-jets, depending if x° — * Xi^b 
occurs. By searching for collider events with > 4 jets plus large E™ ms , along with > 2 
or 3 6-jets which are tagged through the micro- vertex detector, SM backgrounds can be 
reduced by large factors, at only a small cost to signal. 
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This may allow Tevatron experiments to search for gluinos with mass in excess of 400 GeV. 
Such gluino masses are far beyond currently published bounds, and would allow exploration of 
a huge swath of parameter space of Yukawa-unified SUSY models. 

In addition, in the case of the HS model where g — > 66x2 ari a large rate, followed by X2 ~~ ^ 
Xi£ + £~ (typically at ~ 3% branching ratio for each of £ = e or /j,), there may be a corroborating 
signal at much lower rates in the multi-6-jet+i?™ lss + £ + £~ mode, where m(£ + £~) < m~o — m^O' 

We note finally that the results presented here in the context of Yukawa-unified models are 
more generally applicable to any model with very heavy scalars, and large enough tan such 
that gluinos dominantly decay via three-body modes into 6-quarks. They are also applicable to 
models with hierarchical soft terms, where first/second generation scalars are extremely heavy, 
and third generation scalars are much lighter; some references for such models are located 
in EDI. 
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